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CALCULATION OP TURBULENT EXPANSION PROCESSES* 


By Walter- Tollmien 
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./.--On! the basis, of certain formulas recently established by 
L. -Prandtl for the turbulent ' int er change of momentum in sta-r *• 
tionary flows (reference l)-, various eases of "free turbu- 
lence" - that is, of flows without boundary walls - are treated 
in the present report. Prandtl puts the apparent shearing ■ 
stress introduced' by the turbulent momentum interchange 


V - pl ' 


du du 
dy dy 


where 


u average velocity in • ‘x di-recti on 
y coordinate at right angle to x 

l mixing length 

% 

The underlying reasoning is as follow's: The fluid bodies en- 

tering right and left through .a fluid layer with the time av- 
erage value of the velocity u, at turbulence, have the aver- 
age velocity u •*■ l or u - l ££, while the transversely 

... - . dy 

directed mixing velocity is l — , discounting a constant of 

...dy I 

proportionality included in the more or less accurately known 
t of formula (l); 1 , is no constant - at. a. wall 1 = 0. 

The previously cited report by Prandtl (reference l) contains 
a lucid foundation for formula (1). 

, • * * * . * . / f ' 

The present report deals first with the mixing of an air 
stream of uniform velocity wit h the adjacent still air, then 

*"Berechnung turbulenter Ausbr eitungsvorgange . " Reprint 
from Zeitschrift fur angewandte Mathematik und Mechanik, vol . 

6, 19,36 , pp. .1-12. 


* •/. : 
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with the expansion or diffusion of an air jet in the surround- 
ing air space. 2rxpe,r ie-nce indicates that, the width of the 
mixing zone increases linearly with x, if x is the dis- 
tance from the point where the mixing starts. This fact is 
taken int o .account, by" -the ; fpr mu-la 

.v : I = cx (2) 

The constant of proportionality c can as yet he determined 
only hy comparison with experience; it is the only empirical 
constant of the theory. In many instances it will he expedi- 
ent to introduce T] = y/x as a second coordinate. 


1. MIXING OP HOMOGENEOUS AIE STREAM WITH THE 
ADJACENT STILL AIR 

(Two-dimensional problem of the free jet boundary) 


By reason of the limiting .conditions for the average ve- 
locity the formula is preferably expressed with 


u = f (y/x) = f (T) ) 
Then the stream' function is . 

^ = 7 f (.-) dy 


hence 




fCn) d.T) = xP(n) 


3 \|/ 


▼ * — = -r(Tl) + T) P 1 (T) ) • 

OX 

Quantity T is put according to formulas (l) and (2) 


(3) 


(4) 


I = c 3 x 3 


du 

dy 


du 

dy 


The following boundary conditions exist: At the first 

boundary TJ x (homogeneous air stream), u = constant or by 
introduction of a suitable scale u = 1; that is, 
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» , '(T) 1 r . i 


(5) 


furthermore 


: r .* i' i 


■ 'Ml 
I’r »• ’ ■' ‘ Sf] 


^..0 


a conditi o» : by -whi'ch the" continuous connection is secured - 
that is, 


I* " ( T| j ) = 6 


( 6 ) 


' V^):** 0 


that is. 


f(T| h) = Tli 


(7) 


at the second boundary Tjg (still air) must be u = 0; 
that is, 

:•> j" (T1-) = 0 ’ ■ *' 


( 8 ) 


and, to assure , cqnt inuous connection ss’O; that is, 

; a-Q • ' • 


F'»(Tl a ) = 0 

:ev >: c * 5 ‘ * J 1 * 

Since the pressure, in first approximation, can he 

assumed to he constant, tbe equation of motion reads 

:: + ;T * .** , *> 

M + v M = I 


(9) 


dx 


8y p 9y 


Counting y and T| from the still toward the moving air, 
gives, after introduction of the formulas* ■, the equation of 
motion: 


FF» + J3c 0 


. • (1C) 


♦It is readily apparent at this point, that the formula 
u = f ( T| ) ., nece ssarily requires an l proportional to x. 
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which 1 8 solved by F w =0' or F + 2c 3 F ^ = 0. It affords 
uniform velocity in the one case variable velocity in the 
other. The latter solution obviously applies between T1 x 
and T| a , the former, outside of these limits. In the bound- 
ary points the solutions coincide with discontinuity in F " . 
In order to determine the velocity distribution in the mixing 
sone, the di f f erent ial • qquat i on of . the third order 


F * 2c 3 F" =0 (11) 


must be solved. For the time being, a new scale for T1 is 
advisable, so that formula (ll) simplifies to 


F + F" 1 = 0 


(11a) 


The result then is 


F = Cj, e 


-T1 


c 2 e 


T)/ s 


cos 


J3 


T1 + C 3 e 


Tl/s 


sin 


S3 


The five boundary conditions define the constants of integra- 
tion c lt c 2 , c 3 , and, in addition, the still unknown bound- 
ary points T) i and Tl a themselves; 

The calculation is suitably arranged as follows: Put 

T) * T1 - T)j_ or - T) ■ T) + T} x 

so that 

F = d x e“^ + d a e^/ 3 cos T| + d 3 e^/ 3 sin 7) • t . 

2 2 

f 

The boundary conditions are: 

P(tF) = Til, y '(W) = 1, F''(n) S 0 for Tt = 0 , F * ( Tf ) = 0 , 

i * 

F''(tD = 0 for T) a (I2a-e ) 
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From (12a to.--.c-), dj , d a , and d 3 . can be linearly expressed 
in Tlj , equati on. (l 2d), yields T) x * expressed by I) a , ' ,an.d . 
(l2e) finally gives a transcendental equation for fj a , solv- 
able by successive approximation’. It follows that 


Tl a = -3.02, T| x •» 0.981, Tig = -2.04, & x » -0.0062, 

d a » 0.987, d 3 » 0.577 : ‘ 

With this F and F* are. defined as function of TJ . For 
comparison with experience, the original scale must be used - 
that is, t he r educed T| , -.employed thus far, must be multi- 
plied by %/ 2c a , and the reduced F* by U, the velocity of 
the homogeneous air stream. -The curves of the velocities and 
of the streamlines are given in the table and in figures 2 to 
4. The streamline s are plotted for equidistant values of the 
stream function. The streamline emanating from x = 0 is a 
straight line with angle of inclination 


-tan -1 (0.19 V2c a ) 


For comparison, a Gottingen measurement (reference 2) of the 
dynamic pressure distribution with an automatic pressure re- 
corder at the edge of the nozzle of the big tunnel was em- 
ployed. The distance from the nozzle edge was 112 centime- 
ters, the dynamic pressure of the undisturbed Jet q » 56 
kilograms per qquare meter. It can be presumed that the as- 
sumptions of the two-dimensional problem hold .good for the 
Sl,ze of the nozzle. In; figure 5 the calculated dynamic pres- 
sure is shown as a dashed line over the measured dynamic 
pressure. The unknown constant of proportionality c follows 
from the conversion factor for T) , which is = 0.0845. 

The. iwidth of the mixing "zone is 

b = $4c a X 3 . 02 x = 0.0845 X 3.02.x = 0.255 x 

■ . * ! ^ > ... ....... < . . *\ . * t ‘ 

giving a mixing length, of ; r ..- 

I = 0.0174 x * 0.068.2 b • . 
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The relative sma;ll'’n5efrfs s*of l is unusual. The agreement be- 
tween the theoretical and the measured average velocity dis- . 
tribution is very gc^od, 

• vc ? ; • 

2. JET EXPANSION AS TWO-DIMENSIONAL PROBLEM 


Visualize a wall with a narrow slot, which for the study 
may be regarded as being linear, through which a jet of air 
is discharged and mixes with the surrounding still air. As- 
suming, for the first, that the pressure in the jet is the 
same as outside, the application of the momentum theorem af- 
fords a ready separation of the variables - that is, x' and 
T). By reason of the constant pressure the momentum in x 
direction must be 


+OD 


/- 
— 00 


dy m constant 


Putting u a cp(x) f(Tl) results in 

+00 

cp a (x) * J* f S (T|) dT| a constant 


Consequently cp(x) = -i- 

>/~x 


* = -7E f(T ^ J 

vvx 

\|/ * /*^= f(T))dy =/x P(T|) 


(13a) 

(13b) 


V a — p(l)) + -L P 1 (T) ) 11 

2<Jk */x 


(13c) 


The equation of motion can be set up again, which now, 
however, can be immediately integrated once. This interme- 
diate integral can equally be obtained direct from the momen- • 
turn theorem. By marking off a control surface conforming to 
figure 6, tho impulse entering through the lower boundary is 
puv, while on the other side the impulse variation 
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dx 


/■ 


t . 


3 dy occurs. The turbulent shearing stress . 

• • r . • ' c- i i • ; 


T . pc 3 X 3 |a 3a. 

. a . . Sy _ dy 


acts as outside force, hence the relation exists 


uv + [ u a dy » 

OX ^ 


. oo 


Pr.oij this follows, the equation for F (T| ) s 

( v ; :• 

2c 2 F" 3 = FF* 


(14) 

• n * <1 , . r 

r. v •- ■: 4 . v - 


(valid for positive T| , reflected velocity distribution for 
negative T) ) . With a suitable . scale for T) the differential 
equation Is simplified to 


ff« 3 = FF 1 • f ' * "(14a) 


The order of ^hi6 differential equation can be lowered by in- 
troducing z = In P; that is,' P « e* as new dependent 
variable. Then, (z" + z ,3 ) a = z' whence, after putting 
zl,.-=» Z, finally follows 'the differential equation of the 

fir atr.. q?dar : t .. ? v . ' ; 



The solution of the original equation then requires only 
squaring and removal of the logarithms. 

r 

X . • • " ' 

The following conditions must he satisfied: ?or T) = 0 
(center of jet), v * 0 - that is, P * «* * 0. Since 

u ~ P * '* z'e* is not to disappear for T| = 0, z' must be 
of the same ord ; er, of * <» f<jfr T> = 0\ as e* .is of 0. . 5.ejice, 
for TJ ’= 0 by suitable scale determination. 


. .. t 



P = 0 


• ! / 


P' = 1 


(15a) 

(15b) 
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Thus there are afforded two conditions through which the z, 
that satisfies an equation of the second order, Is completely 
defined. • f • - : • ••}**.. 

Th.e "boundary point TJj. itself then follows from the 
condition 

u = 0; that is, z % = Z = 0 (16) 


for the "boundary 1\ r . 

Integration of equation (14) gives 

1 ) . 0 - £ fin C/I + 1 ) - In 
3 l 


The constant of Integration 
z • a Z » 00 for T) = 0: 

0 = 0 - #ys !L at C = -IL 
3 3 7r? 


The condition (16) Z » 0 fot T) r ' yields 



(17) 

C follows from the condition 


0 - -|*/3" tan” 1 (- — ^) * i — = 2.412 
3 . V 3 ./s' 


V3’ . 


4 • 1 • S f * 

Compliance* ‘with equation (15) is predicated" on a study of the> 
behavior of equation (14) for T) » 0, Z * ». As the solu- 
tion (17) in this range is inconvenient, a new form of solu- 
tion which applies for T) = .0; Z * oA"is derived. 


Por Z — ><», obviously 



-Z 2 , that is, Z 

*• . 



hence ''*!'* 

z — > ln T) + c*, * so that P* * z* e* = e Cl ; for * T) ■ 0. 
Thus the last constant of integration follows from equation 
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(15b) as cj = 0; and the :asympt otic approximation for 

T) = 0 follow 8 at - 

; - : • * ' J V • 



0.4 Jl\ + 0.01 Tl a . . 


Z 


in I) - 

3 


+ £^n 3 . . 

3 

(18) 


The quality of the asymptotic approximation (18) is easily ap- 
praised by a comparison with the exact solution (17) in a 
tone in which both forms of solution are appropriate. 


The method is as follows: Compute T](Z) and hence 

Z(T)) a z*(7)) by equation (17), thus obtaining s('H), possi- 
bly by graphical Integration, where in the region about T) • 0, 
the previously determined asymptotic approximation (18) iq- 
taken into account and s' a <*>. Then the desired functions 
P a e z and F' a z'e* are Obtained, by removal of logarithms 
and multiplication. 


The solution P a constant joins the just-derived solu- ; 
tion with :a dlbcOfttinuity in P tt * t oward .the out side . In. the. 
center ("H = : 0) ;F * Act.s as 1 - 0.4 T) 3 ' 8 , which .entails a' 
disappearance of ; the radius of curvature. * For comparison 
with experience, it is necessary to revert from the reduced 
to the original quantiti es a s shown in the table. The conver- 
sion "factor for T| is V~2c' 3 ; the letter s, in the table, 
signifies a characteristic distance from the gap, where the 
speed, in the center of the Jet equals U 8 . According tO ( .(l3a) 
the J 8^eed ' at jet center distant x from the .gap is then 


;., m > v K ( x >.v,yf. 

(See table.) 


3. JET EXPANSION AS R0TA.TIONALLY SYMMETRICAL PROBLEM 


The corresponding rotationally symmetrical problem, in 
which a jet of air discharges, from a very narrow hole in a 
wall, is treated in exactly the lame manner as the two- 


♦ Prandtl has given a refinement of the theory by which 
t-he'' di s appearance of the radius' of curvature in the c'ent.er 
can be avoided. But, since it would lead too far' af ieid , it 
is not discussed here. 


10 


NAOA TM No . 1085. - ?L 


dimensional problem. JVtT»fcV't1ie variables- x :and T) are 
easily separated again. For, on assuming that the pressure 
in the Jet is constant, 

, +oo 


whence for u 


r f 

■I 


Putting 


affords 


tt J* u 3 y dy 
— 00 


y = constant 


. . -i • i *■' t" ■' .• 

u- * I t(n) 

•» r X 


/ 


f(n) ti d-T) * F(n) 


F 1 . F • F 

u = — , v = — - — 

xl) x xTJ 


The differential equation for F is again obtained by inte- 
gration of the equation of motion or by a second Application 
of the Impulse theorem in analogy to figure 6: 


c s (;F" -^p) a . ! ; 


(19) 


Vith the .introduction of a suitable scale for T), the 
differential equation is simplified to 


^F" - II = Ft 1 " 


(19a) 


By substitution: 


there is afforded 


In F, F * e ! 


+ z ' 2 - Ji ) 8 = s' 


and lastly, after Introducing Z = z', the differential equa- 
tion of the first order.. 
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z x * i - z a - yr (20) 

Tl . - ’ .. 

In addition, the. following conditions hold for T) ■ 0: 
u- may not di sappear; . while v. » 0; that is, .'F(o) = e*( 0 '=0, 

while — = z 1 — remains finite and becomes equal to unity 

T] : . ■ .<’• - v j v i >’ ■ . s ii 

by appropriate regularization. , . ? ? J 1 J 

Now a series development of Z(T)) for T| = 0 can be 
applied in such a way that- these conditions are satisfied; z 

must be negative « for T) m 0, in order that e* = 0, 
which is like In T) s ; because P*/T{ then assumes precisely 
a finite value.- ‘ The result is the following development in 
powers of Tl 3 / a t ■ 

2 « | + a A + M| s + cT) 7 / a + d7] 5 + eT| 13 / 3 . . . (21) 


The coefficients are obtained by introduction of th.i-s f-ormii-la/. 
in the differential equation and* comp-ari son^of equal* powers*. 


« = - |V2-. 



JL, 

1715 


37 

240100 * 


e = 0.000014 . . . 


The convergence is poor on approaching the boundary point 
T) r (ZsO), but a development particularly suitable near T) r is 
as follows: Put 

W * Ti r - r\ 

and x •• r . : . * 

Z a a ii 2 + M 3 + c' Tj" 4 + d Tl 5 + "e iT 6 + f Tj" 7 . . . 

. i .■ >• 

and obtain 

a = 1, b = - — t = - d = -i 

4 8 Tjr 64 Tl r a 64 128 J \ T 3 

T _ 19 133 T _ 0.00278 . ... 

e ~ . .i — — ■■■■■ — ■■ IS— ■ ■- *r '■ — • . • 

256.X 5 T) r . 256 X 40T| r 4, T) r a •• 
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The' unknown constant of Integration- T| r is obtained by making 
the values for Z, as known. from the two developments, agree 
in a certain junction point. It results in = 3.4. 

Quantity ' • F • /TV acts • like 1 •- .0.202 T\ 3 ' 3 - in the center; 
the outward junction in F again takes place with a discon- 
tinuity in ' F The conversion factor from the . r/educed to the. 
actual quantity T) is ^/c a ; 8 signifies in the- table a char- 

acteristic distance from the discharge hole for: which the 
speed in the center of the jet is U 8 . 

The computed velocities were compared with Gottingen test 
data.:- (See reference 2.) The diameter of the discharge noz- 
zle was 137 millimeters. The velocity distributions at 100.; 
centimeters and 150 centimeters distance from the nozzle edge .. 
were used for the comparison. This nozzle distance/ a . may *•; 
not be put equal to x, in view of the point discharge ori- 
fice assumed in the present calculation; x is rather com- 
puted from* a by addition of a constant quantity e ' which 
results - for example - from the fact that for greater a, 
for which the. compari son with these calculations is solely 
permissible, the central velocity decreases as 1/x. In the 
present case e = 26 centimeters. Figure 13 shows the theoret- 
ical and the experimental dynamic pressure for a = 100 centi- 
meters, it amounts to 104 kilograms per square meter at the 
discharge orifice; the agreement; of the average values is good, 
aside from a certain asymmetry of the jet which must have had 
different reasons. From the conversion factor for T) follows 

= 0.063 


The radius r of the jet is 

r a 3 / o 3 3.4 x = 0.063 x 3.4 x = 0.214 x 

The mixing distance l is « ex = 0.0158 x = 0.0729 r. 

Zimm (reference 3) has made corresponding experimental 
investigations at considerably lower speed. His findings 
would yield 


= 0.080 


with a dynamic pressure of 5.1 kilograms per square meter in 
the discharge orifice. According to it, a slight increase in 
mixing path by decreasing Reynolds number is likely. 
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4. PREDICTION OF PRESSURE DIFFERENCE 


So far, all cases had been premised on constant pressure. 
This first approximation can be improved by analysis of the 
pressure differences due to impulse variation on the basis of 
the computed speeds and stresses. For the first step, start, 
say, with the second equation of motion, which in the first 
two cases- reads 


' and / 


u ?r + v ~r * — (JM- + 

dx Oy p \ 3 x a y / P Oy 


. i 

3 v / p 


u |r + T |i.iC|i + i d( °yy ) . ft) . i |e 

3x oy • p \ox t y oy y / P oy 

in the rotationally symmetrical case; a y and cr^ are nor- 


mal stresses, respectively effective in y direction or at 
right angles to y and x. Then, integrate with respect to 

Mo * ^/ UTdy - h hf Ti7 - 1 HI * - p HI- ■ - (82) 

0 o 

and 

hi fy**-* i Hi 

0 0 _ „o 


- i Pst dy4 i /!* 

pj y p J y 


d y * - — | 
p 


(23) 


which is equivalent to applying the impulse theorem. If, as 


heretofore, the normal stress, in this case 


and at » 


are discounted, there is obtained 





!>i 


f*2 F F 1 X\ - F 3 ] 711 - 2 / 

^F • 3 T\ d T1 * i [pi 

L K «/ 

P L J 


for the free jet boundary. 
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for axially symmetrical , jet- expansion. With p r denoting the 
pressure at the jet bbundary," p m the pressure at jet center 
and of the homogeneous air stream, respectively, particulari- 
zation of the above formulas yields 

— ~ Pr '= 0.410(2c a ) a/,3 TT a and ? m - ~ - Pr « 0. 248( 2c a ) a/ V a ( x) 

P p 

and 

--y *: •o.ai6(c a ) a/ ® u m a (x) 
p 

Quantity U indicates the speed of the homogeneous air stream, 
and U m (x) the central speed at x. In the first and third 
case, c has been determined, giving 

P m - P r * 0.00584 £JL! 

1 tt \ 

Pm - Pr * -0.0025 £ — flL.?' 

2 

It is apparent- that the thus computed pressure differences, 
being small, .do not cause a substantial modification of . the 
velocities. 

When computing the pressure difference with respect to 
still air, it 'should be borne in mind that at the jet boundary 
a negative pressure equal to the dynamic pressure of the ra-. 
dial inflow speed prevails. With p 0 as the pressure in still 
air 
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*m - P 0 


Pm “ Po s 


0. v 338(.2c S ) 8 ^ 3 p U 3 = 0.00482 1 

' ' ; 2 

0.124(2c*) 8//3 p U n 3 (x). 


: Pm - Po = -0.372(.c a ) a ^ 3 p U m a (x) = -0.00295; £--- m - (x) 

Hence there is positive pressure within the jet in the two- 
dimensional cases, hut negative pressure in the axially sym- 
metrical case. Thi 8 surprising re suit, which also is at vari- 
ance with a rough impulse consideration, points to a defect in 
the theory. The necessary extension will he given in the fol- 
lowing. 


5. EXTENDED THEOREM FOR THE APPARENT STRESSES 


The theorem applied up to now to the stresses introduced 
hy the turbule.ajb impulse exchange 

T = l 3 I |2i •§£, cr = CT= a t = 0 

j dy dy x 7 . , 

• ; . ’* ‘ \ ; • J 

i 8 no more than a first approximation. In any case, it .. pjan. ; he 
easily proved that du/Sy in the cases in point is great with 

respect to . -3X and ^X, - hence the theorem for the mix- 

dx dx dy 

ing speed caused hy the speed difference is good. 

1 1 dy •• ... 

Sp. ih a paturaj, general! zatipn v . of the pr-evious.. theorem, the 
stress tensor is" put equal to ' 


s 


(72 + ▼ V ) 


( VZ • affinor of v; vv is thi'i conjugate affinoP.) r 

'■This relation is important for the calibration of pitot 
tubes in a jet discharging from a nozzle. 


* 
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The stresses to he newly added here, are, in general, neglected. 


except C-, 


au 

j — and 

a* = 2 I s 


by 

ia y 

b 

dy 


X which are used 

y 


for calculating the pressure differences. This port i <5h r.can cels 
in the model problem worked out for the two-dimensional case 
because of the employed boundaries, but not for the axially 
symmetrical case. Here the pressure differences are augmented 
by the integral ’ 

7 


I 


y 


d y 


so that 


Pm ~ Pr 

and 

Pm “ Pq 


+0.15l(c a ) S//3 p U 2 (x) = +0.0012 P 

2 

0.095(c 3 ) 2//3 p U m S (x) = 0.00075 

2 


that i 8 , positive pressure within the jet, as in the other 
cases . 


Translation by J. Vanier, 
National Advisory Committee, 
for Aeronautics. 
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Figure 13.- Calculated and recorded dynamic pressures for a ■ 100 
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